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Abstract 



o 

We study a quaternion version of determinantal random point fields denned 
by using the Moore-Dyson quaternion determinant. First, we give sufficient 
conditions that ensure that a self-dual quaternion kernel defines a valid ran- 
dom point field, and then we prove a CLT for quaternion determinantal fields. 
The proofs are based on a quaternion extension of the Cauchy-Binet determi- 
nantal identity. 

Oh 

1. Introduction 

A determinantal random point field is a random collection of points on a nice topological 
space with the probability distribution that can be written as a determinant. The determinantal 
' point fields describe various mathematical objects including eigenvalues of random matrices, 

zeros of random analytic functions, non-intersecting random paths, and spanning trees on net- 
works. It is conjectured that they are also related to other important objects such as Riemann's 
zeta function zeros and the spectrum of chaotic dynamical systems. See E4l . |fT3ll , Ifl2l and 
101 for reviews. 

The definition of the determinantal point fields uses the standard matrix determinant. How- 
ever, in some applications, the distribution of a random point field can be represented as a 
determinant of a quaternion matrix. An important example is provided by eigenvalues of or- 
thogonal and symplectic random matrix ensembles. It is natural to call these fields as quater- 
nion determinantal fields and study their properties. This paper makes some initial steps in 
this direction. 

In order to define a quaternion determinantal point field, we need to establish some pre- 
liminary notations. 

A random point field X = (X, B, P) on the set Acl" (or Z n ) is a probability measure P 
on the space X of all possible countable configurations of points in A. It is convenient to think 
about this object as a collection of functions that sends every n-tuple of non-negative integers 
(ki, . . . , k n ) and every n-tuple of Borel subsets of A, {A\, . . . , A n ) to a non-negative number, 
which can be inteipreted as a probability to find ki points in the set Ai . These functions are to 
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satisfy some consistency conditions, which we do not specify here. The reader is advised to 
consult paper [14] or book [6] for more detail. 

Let # (A{) denote the number of points of X located in the set A^. 

Definition 1.1. A locally integrable function R k : A fc — > is called a Appoint correlation 
function of a random point field X = (X, 23, P) with respect to background measure /J,, if for 
any disjoint Borel subsets A\ , . . . , A m and any non-negative integers k\ , . . . , k m , such that 
YliLi ki = k, the following formula holds: 

m ,. 

EjJ[#(i4 < )...(#(i4 i )-fcj + l)]= / R k (xi,...,x k )dn(x 1 )---dn(x k ), 

(i) 

where E denote expectation with respect to measure P. 

Note that on the left is the expected number of ordered configurations of points such that 
set Ai contains ki points. Note also that the correlation functions are defined only up to sets 
of measure 0. 

Definition 1.2. A random point field X on the set A C W 1 will be called a quaternion deter- 
minantal point field if its correlation functions can be written as quaternion determinants: 

R m (x 1 ,... ,x m ) = DetM (K (ari , Xj- )) | i<i, 3 -< m ' ( 2 ) 

where K (x, y) is a self-dual quaternion kernel (that is, K (y, x) = (K (x, y))*), and Det-M 
is the Moore-Dyson quaternion determinant. 

(We will give more information about quaternions and quaternion determinants in the next 
section.) 

Example 1. It was shown by Dyson [8 ] that the eigenvalues of the circular orthogonal and 
symplectic ensembles of random matrices form a quaternion determinantal field. Later this 
result was extended by Mehta (see Chapters 7 and 8 in [ 18']) to the case of Gaussian orthogonal 
and symplectic ensembles. We will say more about the circular symplectic ensemble in the 
last section. 

Example 2. (Ginibre on Q) Let z and w be in Q, where Q denotes quaternions over the 
real field, and define the kernel 

*«<-> = £^ 0) 

This kernel defines a random field on the set of all real quaternions Q ~ M 4 taken with the 
background measure dp, (z) = 7r~ 2 e - ' 2 ' dm (z) , where dm (z) is the Lebesgue measure on 
Q. The fact that kernel ((3]) defines a random point field follows from Proposition @3]below. 

Recall for comparison that the usual Ginibre random point field on C gives a distribution 
of eigenvalues of a random n-by-ra complex Gaussian matrix. It is determinantal with the 
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kernel 

K.M = ± 

k=0 

where z and w are in C and the background measure is dfi (z) = Ti- l e~\ z \ dm (z) . (See ifTTH . 
Section 15.1 in |QJD, and Section 4.3.7 in [12).) 

Example 3. (quaternion Ginibre on C) This is another generalization of the Ginibre 
random point field. See Section 15.2 in |fl8l for details. Let 

^ v ) = J_ y 111 1 ( n _ 

ri v ' ; 2tt ^ 2H\ (2j + 1 v ; 



and define the quaternion kernel by its complex matrix representation: 

<p(K N (z,w)) = 



(f) N (w, z) 4> N (w, z) 
(j> N (z,w) (p N {z,w) 



(The definition of the bijective map tp: Qc — > M2 (C) is standard and given in the next sec- 
tion.) Then this kernel defines the determinantal kernel with respect to the signed background 
measure d[i (z) = e - ' 2 ' (z — z) dm (z) . 

Example 4 (Bergman kernel on Q). Let z and w be in Q, and define the kernel 

n 

K n (z,w) = ^2(k + 2)z k (w*) k - (4) 

fc=0 

This kernel defines a random field on the unit disc of real quaternions with the background 
measure dfi (z) = ir~ 2 dm (z) , where dm (z) is the Lebesgue measure. 
For comparison, the Bergman kernel on the unit disc in C is given by 



(w) k 



K n (z,w) = y^{k + l)z k 

k=0 

and corresponds to the distribution of zeros of power series with i.i.d complex Gaussian coef- 
ficients. (See Section 15.2 in (Hi). 

The quaternion determinantal point fields are a particular case of the Pfaffian point fields 
which were studied in ED, E5l and 131 . 

Which self-dual quaternion kernels define random point fields? 

Note that a correlation function is automatically symmetric if it is defined as in d2]). That 
is, we have 

(-£(t(1) j • • • > ^ir(m)) = { x l 1 • • • > -^m) 

for every permutation a G S m . Hence, by the Lenard criterion ( |[T5l and |[T6l ). positivity is 
a necessary and sufficient condition that ensures that a kernel corresponds to a random point 
field. This condition can be explained as follows. Let X be the space of configurations of 
points in A. Let cp = {tfk} denote an arbitrary sequence of real-valued functions tp^ over 
A k which have the property that they are zero if at least one of its arguments is outside of a 
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compact set K in A. Define operator S on tp as follows: S maps sequence (p to a function 
onl 

where x = (xi,x%, . . .) is an enumeration of a point configuration in X, and the sum is ex- 
tended over all finite sequences of distinct positive integers (i±, . . . including the empty 
sequence. (This sum is essentially finite since the definition of the space of configurations 
requires that the number of points of any configuration in any bounded subset of A be fi- 
nite.) The positivity condition says that if (Sip) (x) > for all x E X (including the empty 
sequence), then it must be true that 

/ p(x)dp := ip + Y] / (xi,.. ■ ,x k )R k (xi, . . .,x k )dp(x 1 ) . ..dp(x k ) > 0. 

However, this criterion is often difficult to verify. 

In the case of the usual determinantal fields with self-adjoint kernels a much simpler crite- 
rion is given by the Macchi-Soshnikov theorem (see ifTTl and |[24l ) that says that a self-adjoint 
kernel K (x,y) defines a determinantal point field if and only if the corresponding operator 
K is in the trace class and all its eigenvalues are in the interval [0,1]. This is equivalent to the 
requirement that both K and I — K we. non-negative definite. 

Our first result is a weaker version of this theorem for quaternion determinantal fields. 

Suppose that a quaternion kernel K can be written as follows: 

oo 

K ( x , V) = X kU k (x) u% (y) , (5) 

k=l 

where are scalar and u k (x) is an orthonormal system of (possibly complex) quaternion 
functions: 

u* k (x) ui (x) dn(x) = 5 ki . 

A 

(The series in (f5]) are assumed to be absolutely convergent almost everywhere.) We will say 
in this case that K (x, y) has a diagonal form with eigenvalues Afc. If X k are real and u k (x) is 
an orthonormal system of real quaternion functions, then we will say that K (x, y) has a real 
diagonal form. 

Theorem 1.3. Suppose that a quaternion kernel K (x, y) has a real diagonal form with eigen- 
values X k . Assume that all X k G [0, 1] , and that Y2k=l ^ k < 00 • Then the kernel K (x, y) 
defines a determinantal point field with finite expected number of points. 

The conditions of this theorem are sufficient but not necessary since there exist self-dual 
quaternion kernels which do not have a real diagonal form and still define a valid point field. 
For example, let space A consist of two points and has the counting background measure. We 
can define a quaternion determinant point field on this space by setting the kernel equal to a 
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2-by-2 quaternion matrix 

K = 

where a = (3ii— 5j) /4 so that a 2 = —1. This matrix is the matrix of the orthogonal projection 
on vector (1, a) . Clearly, this matrix is self-dual with determinant zero. Hence, it defines a 
random field that has exactly 1 point uniformly distributed on A. On the other hand, this kernel 
has a diagonal form with real eigenvalues but it does not have a real diagonal form. 

Moreover, there exist self-dual quaternion kernels which do not have any diagonal form 
at all and still define a valid point field. Here is an example of this situtation. Again, let the 
space A consist of two points, and consider the kernel equal to a 2-by-2 quaternion matrix 



K 



1 a 
-a 1 




where a = ii — j, so that a 2 = 0. Clearly, this matrix is self-dual with determinant 1. Hence, 
it defines a random field with R\ = 1 and the correlation matrix 

R2 = 

This random field has exactly 2 points. On the other hand, K does not have a diagonal form 
because if it had such a form then it would be a projection matrix, but an explicit calculation 
shows that K 2 / K. 

(Note by the way that the determinantal field defined by K is the same as the field defined 
by the identity matrix. Hence, two different quaternion kernels can correspond to the same 
random point field.) 

On the other hand, in Theorem 11.31 it is not possible to omit the requirement that the 
eigenfunctions (x) are real quaternion. Indeed, there are kernels that have a diagonal form 
with real eigenvalues in the interval [0, 1] but that do not define a random point field. Condider 
for example, the same two-point space A and the kernel 




It has a diagonal form with A = 1 and u = (2/v3j [1, i/2}* but it does not define a random 
field since the first correlation function is negative at the second point. 

Most of the examples from random matrix theory are concerned with the kernels that do 
not have a real diagonal form. Orthogonal random matrix ensembles correspond to kernels 
without diagonal form and symplectic ensembles correspond to kernels without real diagonal 
form. Hence, it appears desirable to improve Theorem 1.3. What are sufficient and necessary 
conditions for a self-dual quaternion kernel to define a determinantal random point field? 

While this question is not answered in this paper, we can slightly generalize Theorem 11.31 
to give sufficient conditions for kernels without a real diagonal form. We will say a kernel 
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K (x, y) has a quasi-real diagonal form (or simply call it quasi-real) if it has a diagonal form 
with real eigenvalues. We will call it positive if 

Det M (if (a*,a:j))li<ij<m > 

for all m G Z + and all x\, . . . ,x m . It is remarkable that quasi-real kernels with positive 
eigenvalues are not necessarily positive (as the example above shows). 

We will call a quasi-real kernel K (x,y) completely positive if every of the kernels 

Ki = u i ( x ) u *i (v) 

is positive, where {ui (x)} is an orthonormal set of eigenfunctions of K(x,y) and I = 
[i\ , . . . , i m ) denote an ordered subset of indices of all eigenfunctions. 

Theorem 1.4. Suppose that a quaternion kernel K (x, y) is finite-rank, quasi-real and com- 
pletely positive. Assume that all A& 6 [0, 1] . Then the kernel K (x, y) defines a determinantal 
point field. 

Remark: We omit the assumption that J2h=i < oo which we imposed in Theorem II .31 
since we assume that the kernel is finite-rank. It should be possible to extend this theorem to 
a more general case when there are infinite number of A& and YlkLi < o°- 

Still, even if we restrict attention to kernels with a diagonal form, the conditions of The- 
orem 11.41 are not necessary. For example, consider the two-point space A with the counting 
measure p,. Let 

s /19 20 \ . 
a = 1 + 2i + i i, 

and define the kernel by the matrix 

K = 

Then Det^/ (K) « 0.3745 and therefore the kernel defines a valid random point field on A. 
On the other hand the eigenvalues are Ai,2 ~ \ ± 0.3529i and therefore the kernel is not 
quasi-real. 

More generally, we have the following example. 

Let A be a complex number and assume that | A| < 1, and ReA > |A| 2 . Let u (x) and v (x) 
be two orthonormal quaternion functions on A and let K (x, y) = u (x) u* (y) + v (x) v* (y) . 
Assume that that u (x) u* (x) = v (x) v* (x) for every x and that it is true that at 

/ K{x,x) K(x,y) \ 

for all x, y G A. 
Then the kernel 

An (x) u* (y) + Xv (x) v* (y) 
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corresponds to a quaternion determinantal random field on A. Indeed, it is a mix of the 



Numerical data suggest that in the case of the symplectic random matrix ensembles the 
kernel is quasireal. 

The proof of Theorem [T3]is based on the following result. 

Suppose K (x, y) has a real diagonal form with eigenvalues A&, that all A& G [0, 1] , and 
that J2T=i Afc < oo. Define a random kernel 



k=l 

where are independent Bernoulli random variables. The random variable takes value 1 
with probability 

We will prove later that this kernel defines a quaternion determinantal field. 

Theorem 1.5. Suppose K (x, y) has a real diagonal form with eigenvalues A&, that all A& G 
[0, 1] , and that Y2k=i < °°- Let X^ be the random point process which is a mix of the 
quaternionic determinantal processes with random kernels (x, y) defined in (0. Then, Xg 
is the quaternion determinantal process with the kernel K (x, y). 

This theorem is a quaternion version of Theorem 7 in lfl3ll . In order to prove Theorem 1 1.5 1 
we need an analogue of the Cauchy-Binet formula for quaternion determinants which we state 
in Section [3] and prove in Appendix. 

An analoguous result holds for the quasi-real case. 

Theorem 1.6. Suppose K (x, y) is finite-rank, quasi-real and completely positive with eigen- 
values Xk G [0, 1] . Let X^ be the random point process which is a mix of the quaternionic 
determinantal processes with random kernels (x, y) defined in (|6]). Then, X^ is the quater- 
nion determinantal process with the kernel K (x, y). 

Another consequence of Theorem ll.5l is the central limit theorem for the number of points 
in quaternion determinantal fields. 

Theorem 1.7. Let X n be a sequence of quaternion determinantal point fields defined on a set 
DC1 with kernels that have a real diagonal form with eigenvalues in the interval [0, 1]. Let 
M n denote the number of points of X n in D. Suppose that every M n is finite with probability 
1, and thatYar (Af n ) — > oo as n — > oo. Then, the sequence of random variables 



approaches a standard Gaussian random variable in distribution. 

This is a quaternion analog of a theorem that was proved by Costin and Lebowitz in Q 
and Diaconis and Evans in [7] for particular cases, and by Soshnikov in ll26l for general 




oo 




(6) 



M n -E(Af n ) 
VVar (Af n ) 
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determinantal ensembles. Later, a simplified proof was suggested in |[T3l and we use its main 
idea to prove our theorem. 

For the quasi-real case we have the following parallel result. 

Theorem 1.8. Let X n be a sequence of quaternion determinantal point fields defined on a 
set D C R with kernels that finite -rank, quasi-real and completely positive, and which have 
eigenvalues in the interval [0, 1]. Let M n denote the number of points of X n in D. Suppose 
that every M n is finite with probability 1, and that Var (Af n ) — > oo as n — > oo. Then, the 
sequence of random variables 

Af n -E (M n ) 

approaches a standard Gaussian random variable in distribution. 

We will see in the final section that the kernel of the circular symplectic ensemble of 
random matrices is quasi-real, which gives an independent proof that it defines a determinantal 
point field on the unit circle. Numerical evaluations suggest that this kernel remains quasi- 
real even if it is restricted to an interval of the unit circle. However, the proof of this claim is 
elusive. 

The rest of the paper is organized as follows. Section [2] provides some background on 
quaternion matrices and determinants. Section [3] formulates the Cauchy-Binet identity. Sec- 
tion Improves Theorem 1 1.31 Section Improves Theorem 1 1.7 1 and Section [6]provides an illustra- 
tion of the general theorems by using the circular symplectic ensemble of random matrices. 
The appendix contains a proof of the quaternion Cauchy-Binet identity. 



2. Quaternion matrices and determinants 

Recall that the algebra of real quaternions Q is a non-commutative division algebra isomor- 
phic to a subalgebra of the algebra of all 2-by-2 complex matrices (C) which is generated 
over R by the identity matrix and matrices 




which are called quaternion units. The quaternions are typically written as q = s+xi+yj+zk, 
where s, x, y, z are real. The conjugate of q is q* = s — xi — yj — zk, and the norm of q is 

|9| := (g*«) 1/2 = (a 2 + x 2 + y 2 + « 2 ) 1/2 . 

It is also useful to define complex quaternions as linear combinations of the identity el- 
ement and quaternion units i, j, and k with complex coefficients. The algebra of complex 
quaternions Qc is isomorphic to M2 (C). In the case of complex quaternions, we use the 
same definintion for the conjugate, q* = s — xi — yj — zk , and we define the norm as 

/ o i2 i2 ii2\ 

IMI := (I S I +Fl + M +M ) . We say real quaternions for quaternions with real co- 
efficients. 
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In terms of 2-by-2 matrices, if 



Q = 



( 



a 



c 



b 
d 



) 



then its conjugate is 




— c 



d 



-b 



a 



) 



Quaternion matrices are matrices whose entries are quaternions. They obey the usual 
rules of matrix addition and multiplication. The dual of a quaternion matrix X is defined 
as a matrix X*, for which (X*) lk = (X^)* . Self-dual quaternion matrices are defined by 
the property that X* = X. Unitary quaternion matrices are matrices with the property that 
X*X = XX* = I, where I is the identity matrix. 

If we represent each quaternion by a 2-by-2 complex matrix, then matrix X becomes a 2n- 
by-2n complex matrix which we denote ip (X). Let J be a 2n-by-2n block-diagonal matrix 
with the blocks 



Complex matrices A that satisfy the condition A T = —JAJ are called skew-Hamiltonian, 
and the matrices V that satisfy the condition V T JV = J are called symplectic. Hence, a 
quaternion matrix X is self-dual if and only if (p (X) is skew-Hamiltonian, and X is unitary 
if and only if ip (X) is symplectic. 

Some of more involved linear algebra concepts are known to be valid for quaternion ma- 
trices with real quaternion entries, for example, matrix eigenvalues and singular value decom- 
position. A number A is called an eigenvalue of quaternion matrix X if for some non-zero 
quaternion vector v, we have Xv = vX. (These are the right eigenvalues of X, which are the 
most convenient in applications.) It is easy to see that if A is an eigenvalue, then q~ 1 \q is 
also an eigenvalue for any quaternion q. However, for self-dual quaternion matrices with real 
quaternion entries, all eigenvalues are real and it is possible to show that any n-by-n matrix 
X of this type has exactly n eigenvalues (counting with multiplicities); see Zhang ll28l . 

Furthermore, it is possible and useful to generalize the concept of determinant to quater- 
nionic matrices. There are several sensible ways to do this and in this paper we will only use 
the Moore-Dyson and Study determinants. Interested reader can find details in a review paper 
by Aslaksen (H. 

The Study determinant of X is defined as the usual determinant of the complex represen- 
tation of X, 




on the main diagonal. Then it is easy to check that 



<p(X*) 



j ip (xy j. 



Det s (X) := det (if (X)) . 
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It can also be defined in a slightly different way for real quaternionic matrices. If 

X = X 1 + X 2 \ + X 3 j + X 4 k, 

where X\ , X2 , X3 , and X4 are real, then we define two complex matrices A = X\ + X 2 i and 
B = X3 + X^i. Then, tp (X) is a 2n-by-2n complex matrix defined by the following rule: 




1> (X) = _ _ , (7) 



where ^4 is the conjugate of matrix A, that is, (A) kl = (A^)* , and similarly for B. The 
matrix ip (X) is called the complex adjoint of matrix X. Then 

Det s (X) = det (V (X)) , 

where the determinant on the right is the usual determinant of complex matrices. 

The Study determinant is multiplicative: Det^ (AB) = Detg (A) Dets (B) for square 
matrices A and B. 

The Moore-Dyson determinant of a self-dual matrix X with real quaternion entries can be 
defined as the product of the right eigenvalues of the matrix. Remarkably, this determinant 
can also be extended to all matrices by using a variant of the Cayley combinatorial formula 
for the determinant. Namely, let S n be the group of permutations of the set { 1, . . . , n} . Write 
every permutation a as a product of cycles: 

a = (nii 2 ...i s ) (n 2 j 2 ■ ■ ■ jt) ■ ■ ■ (n r k 2 ■■■h), 
where ni are the largest elements of each cycle and n\ > n 2 > . . . > n r . Then we can write 

Det M pO = £e(<7) {X n ii2^"*2«3 • • • -^Q«ni) • • • {X rlj .f :2 Xf (:2 ) c . i . . . Xk irlr ) , 

where e (a) = (— l) n ~ r is the sign of the permutation a. (see [19] and [HI). 

Note that this definition allows one to calclulate the quantity Det-M (X) for an arbitrary 
quaternionic matrix. Dyson established that this quantity is scalar for a self-dual matrix, that 
is, in this case, the i, j, and k components of the determinant are zero. 

The Moore-Dyson quaternion determinant of a self-dual quaternion matrix can also be 
written as the Pfaffian of a related complex matrix. If X is a self-dual quaternion matrix, 
then — Jip (X) is antisymmetric (that is, [— Jip (X)] T = Jtp (X)), and we can compute the 
Pfaffian of this matrix. We have 

Det M (X)=Pf(-J<p(X)). 

(see H ] and Proposition 6.1.5 on p. 238 in Forrester's book [10]). 

In terms of the transformation ip, this can be written as follows. Let X be a real quaternion 
matrix and let A and B be defined as in ©. Let 
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Then, 

-J$ (X) = 

If the real quaternion matrix X is self-dual, then — Jtjj (X) is antisymmetric, and it can be 
shown that 

Det M (X) = -Pf (*)) • (8) 

The Study and Moore-Dyson determinants are related by the following formula: 

Dets(X) = Bet M (X*X). (9) 
(See formula (6.13) on page 239 in [10] or Corollary 5.1.3 on p. 75 in QJQ .) 

3. Cauchy-Binet formula 

Determinants of matrices with quaternion entries have been studied for almost a hundred 
years. (Study wrote a paper [27 ] about them in 1920, and Moore presented his definition |fl9l 
in 1922). In the 1970s, this subject got a boost after Dyson re-discovered Moore's determinant 
and related it to the distribution of eigenvalues of random matrices, (H. In recent years, the 
entire subject of quaternion linear algebra received an increased share of researchers' attention 
because of its applications in image processing ( l20l . |9l , l22l . ifTI). 

Unfortunately, while there are several variants of quaternion determinants, none of them 
enjoys all the properties of the usual determinant. The validity of each standard determinantal 
identity has to be checked individually. 

Recall that the Cauchy-Binet identity states that if A is an m-by-n matrix and B is an 
n-by-m matrix, with n> m, then 

det (AB) = det {A 1 ) det (B 1 ) , 
i 

where the summation is over I = (ii < i% < . . . < i m ), the ordered subsets of {1, . . . , n} that 
consist of m elements. Matrices A 1 and B 1 are square matrices that consist of m columns of 
A and m rows of B, respectively, with indices in /. 

An implicit assumption in this result is that the entries of the matrices A and B are from a 
commutative ring, for example from a field of complex numbers. Unfortunately, in this form 
the Cauchy-Binet identity fails for the quaternion determinants. 

It turns out that a weaker form of the Cauchy-Binet identity still holds for quaternion 
determinants. 

Theorem 3.1. Let C be an n-by-m matrix with (possibly complex) quaternion entries, n> m, 
and let C* be the dual of C. Then 

DetM (C*C) =J2 Det M ((C 1 )* C 1 ) , 
i 

where the summation is over I = (ii < < . . . < i m ), the >ordered subsets of{l, . . . , n} 
that consist of m elements, and where C 1 consists of rows . . ,i m ofC. 
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A proof of this theorem is in Appendix. 

Corollary 3.2. Suppose that A is an n-by-n diagonal matrix with scalar entries Aj on the 
main diagonal and that C is an n-by-m quaternion matrix, n>m. Then, we have. 

Det M (C*KC) = Y, A 'i ' ' ' X imVet M ((C J )* C 1 ) . 

I=(il,...,im) 

i\<...<im 

Proof of Corollary |33| Let A 1 / 2 be an n-by-n diagonal matrix with scalar entries such 
that (A 1 / 2 ) = A. For / = . . . ,i m ) , let (A 1 / 2 ) denote an m-by-m matrix which is 
formed by taking entries at the intersection of rows and columns i\, . . . , i m in matrix A 1 / 2 . 
We write 

Det M (C* AC) = ^DeW ((a 1 / 2 c) J *(a 1 /2 C -)^ 
= EDet,((AV 2 c) 7 ) 
= £Det 5 ((AV 2 )V) 
= ^Det s ((A^)^Det5(^) 
= Yl A il ...A im Det M ((C0*C J ). 

I=(il,...,im) 

The first line is the Cauchy-Binet identity. The second line follows from (O. The fourth 
line follows by multiplicativity of the Study determinant. And in the fifth line we used (© 
again. □ 

4. Quaternion determinantal fields 

Proposition 4.1. Let Kn (x, y) = J2k=i u k ( x ) u t (v) > where (x) are orthonormal quater- 
nionic functions, and assume that Kjy (x, y) is positive, that is, 

Det M (K N (ari,a;j-))li<<j< m > 

for all m G Z + and all x±, . . . , x m . Then Kn (x, y) defines a valid quaternion determinantal 
field. 

Note that for the case when Uk (x) are real quaternion function the assumption of positivity 
can be dropped. Indeed, for every collections of points x±, . . . , x m , the matrix 

K= K r ( Xi ,Xj\ j=l ^ m 

is positive semidefmite. (That is, for every real quaternion vector v, v*Kv > 0.) Hence, all 
eigenvalues of this matrix are real and non-negative. Since for self-dual matrices with real 
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quaternion entries the Moore determinant can be defined as the product of the eigenvalues, we 
conclude that the determinant of K~n (xi , Xj ) is non-negative. 

The proof of the proposition follows the standard lines. We include it for completeness. 

Proof: We need to show that the functions defined by the rule 



R m {xi,...x m ) = DetM (K N (xi,Xj)) , i,j = 1, ... 



in 



are the correlation functions of a random field. 

First, by assumtion of positivity all functions R m {x\, . . . x m ) are non-negative. 
Next, by integrating the kernel we find that 



/ 

Jr 



K]\f (x, x) dx = N, 
K N (x, y) K N (y, x) dydx = N. 



By using formulas (Q]), we can conclude that the total number of points in the process with 
kernel Kn (x, y) is exactly N. (Its expectation is N and its variance is 0.) Hence, it remains 
to show that the functions R m (x\, . . . x m ) agree among themselves for all m. This can be 
done by the quaternion analogue of the Mehta lemma. (Compare Theorem 5.1.4 on p. 75 in 
Mehta flgft 
Let 

K m • — (J^yEii 3?j'j)l<i,j<m 

for a self-dual quaternion kernel K. 
Lemma 4.2. Assume that K satisfies either 

[ K(x,y)K(y,z)dy = K(x,z) (10) 

or 

V U K(x, y)K(y, z)dy^j = tp (K(x, z)) + E<p (K(x, z)) - <p (K(x, z)) E, (1 1) 
where 

*-(;?)■ 

Then 

[ Det M (K m )dx m = (N-m + l)Det M (Jm-i), (13) 
where N = J K(x, x)d/j,(x). 

This result is due to Dyson. (See proof of Theorem 4 in [8].) 
For the kernel Kn (x,y) , equation (fTOl holds, and we find that 



R m (xi, . . .,x m ) dx m = (N - (m- 1)) Rm-1 (x\, ■ ■ .,X m -l) , 
which shows that all correlation functions are all in agreement. 



14 



V. KARGIN 



Corollary 4.3. Assume that 

oo 

K ^(x,y) = ^2^kU k (x)u* k (y) , (14) 

k=l 

where every £ k is an independent Bernoulli random variable that takes value 1 with probabil- 
ity A/%, and (x) are orthonormal real quaternionic functions. Suppose that Y2k=l ^ k ^ °°- 
Then with probability 1, the function (x, y) is a kernel of a determinantal point field, X^. 
For a given realization of variables the total number of points in X^ is Y^k=i Cfc- 

We can prove an analogous result for a quasi-real diagonal form. 

Corollary 4.4. Assume that 

N 

K ^{x,y) = s ^ j i k u k (x)u* k (y) , (15) 
fc=i 

where every £/. is an independent Bernoulli random variable that takes value 1 with probabil- 
ity Afc. Suppose that the kernel K (x, y) = J2k=i u k ( x ) u t (v) JS completely positive. Then 
with probability 1, the function Kg (x, y) is a kernel of a determinantal point field, X^. For a 
given realization of variables the total number of points in X^ is X^fcLi £fc- 

Proof of Theorem ll.5t By the iterated expectation formula, the conelation functions of 
process X^ equal EDet m (K{ {xi, Xj)) , where the expectation is taken over randomness in £. 
Hence, it is enough to prove that 

EDet M (K^(xi,Xj)) =Det M (K (xi,Xj)) , i,j = l,...,m, (16) 

almost everywhere. 
First, let 

R 

Kr (x, y) = ^2 XkUk ( x ) u l (y) ■ 

k=l 

Since we assumed the absolute convergence of the kernel, hence DetM (Kr (xi,Xj)) con- 
verges to DetM (K (xi, xj)) almost everywhere as R — > oo. 
Let i?-by-m matrix C be defined as 

C k i = u* k (xi) , k = l,...,R; 1 = 1, ...,m. 

and let A be an R-by-R diagonal matrix with diagonal entries Aj. 
By Corollary [321 

Bet M (K R (x i ,x j )) = ]T A il ...A iro Det M ((C / )*C / ). (17) 

I=(il,...,im) 

i\<...<i m 

Next, let the random variable Det m (-fQ Xj)) be denoted as Y and let Ar be the event 
that all £k are zero for k > R. (That is, Ar = r\k> R {£,k = 0}). Note that 

EY = E (Y\A R ) P (A R }) + E (Y\A C R ) P (A%) . (18) 
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By using independence of Ar and for k < R, we find that 

E (Y\A R ) P (A r ) = EDet M (C*A € C) P (A R ) , 

where denotes an R-by-R diagonal matrix with diagonal entries £j. 
Next, by Corollary [Ol 

EDetAf (C*AgC) = E £ ^ . . . £ im Det M ((C 7 )* C') . 

I=(il,...,tm) 
ii<...<i m 

Since the variables ^ , . . . , ^ m are independent and have expectation , . . . , Aj m , we 
find that E^ . . . & m = A 4l . . . A im . Hence 

E (T|Ar}) P = Det A f (tfji (a*, xj)) P (A*) , (19) 

and the probability P (A#) converges to 1 as R — > oo by independence of Borel-Cantelli 
lemma and the assumption ^ A& < oo. 

Now let us show that E (Y \ A R }) P (A R ) converges to zero almost everywhere. By posi- 
tivity of the determinant, it is enough to show that 

/ E (Det M (# e (xi, Xj )) \A R ) P (A r ) (20) 

as i? — > oo. 
Let 

OO 

fc=l 

which is finite since both the expectation and the variance of the sum on the right hand side 
are convergent. Since is a projection operator, the total number of points of the process A^ 
in R equals n^. By changing the order of the expectation and the integral signs in (1201) . which 
is possible since the integrand is positive, and by using the identities for correlation functions 
we obtain that we need to estimate 

Efa fa - 1) ... fa - m+ 1) \A R )F(A R ) , 

which is smaller than 

E(nf\A c R )F(A R ). 

By expanding 




and using the fact that = 6c f° r every integer s > 1, we observe that it is enough to show 
that 

e( E ^■■■Ua c r )f(a r ), 

\ii<...<i r / 

where the sum is over all possible ordered r-tuples , , . , i r ) such that %\ < . . . < i r , and 

1 < r < m. 
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We can divide the sum in two parts. The first part is when i r < R. In this case, the variables 
£i 2 , . . . , £,i r are independent from the event A C R , and therefore we can estimate this part of 
the sum as 

]T A n X i2 . . . \ r \A R \ P (A C R ) < cS m W (A R ) , 

i<...<i r / 



\i 



where S := YlkLl an d c is an absolute constant. Hence this part converges to zero as 
R — > oo because P (A C R ) converges to zero. 

The second part of the sum is when i r > R. In this case, the event £j r = 1 implies A R and 
therefore, 

1 Ci T \A c R )¥(A c R ) = F& 1 = l,...,£ ir = l,mdA R ) 

= P(e il = l,... > e< r = l) 
= Ajj . . . Aj r . 

Therefore, we estimate the second part as 



CO 



Aii • • ■ Aj r < Aij . . . Aj r 

tl<...<V, i r =i?+l ii,..i r _i = l 

oo 

i r =_R+l 

Hence the second part converges to zero as R — > oo, because the tail sums Y1T=R+1 A v 
converge to zero. 

This shows that E (y| A^}) P (A C R ) in ( fT8l ) converges to zero as R — > oo. If we compare 
(fTTT ) and ( fl9l ) and let i? grow to infinity, we find that 

EDet A f (K£ ( Xi , Xj )) = Det M (K (x h Xj )) (21) 

almost everwhere, and this completes the proof of the theorem. □. 

Proof of Theorem: 11.61 The proof is essentially the same as for Theorem 1 1.51 but it is 
easier, since we assumed that the number of eigenvalues is finite. The only difference in the 
proof is that the matrix C is defined differently: C is s-by-m matrix defined by 



C k i = ^2 Vkill i ' k = 1 ' • • • ' s; l = h--- 

i=l 

□ 

Proof of Theorem 11.31 

This is an immediate consequence of Theorem 1 1.51 □ 
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5. A CLT FOR QUATERNION DETERMINANTAL FIELDS 

Proof of Theorem [L7j The assumptions ensure that every field X n has a kernel K n (x, y) 
which can be written as Y^k=i ^k^ u k^ ( x ) u t^* (v) w * tn f unct i° ns ( x ) wrnc h are or - 
thonormal on D and are real quaternionic, and with real a£ G [0, 1] that satisfy the condition 
J2T=i ^1 < 00 • ^ * s c l ear th at Theorem 1 1 . 5 1 holds also for quaternion determinantal fields 
defined on D, and therefore we can write 

00 

k=l 

(n) ( (n) \ (n) 

where Q are independent Bernoulli random variables, and P I Q = 1 1 = A^, . We can 
truncate these infinite sums: 

fe=l 

in such a way that M n — M n approaches in probability, lW n — > lW n , YarJ\f n — > YarJ\f n . 
By a suitable modification of the Lindenberg-Feller theorem ([23 ], Theorem III.4.2 on p. 334) 
we have that 

N n - E 

^Yar (AT;) 

approaches the standard Gaussian random variable and this implies the claim of the theorem. 

□ 

The proof of Theorem ll.8l is similar. 

6. An Example 

Recall that the circular unitary ensemble (CUE) is defined as the probability space of all 
A^-by-Af unitary matrices with the Haar measure. The circular orthogonal ensemble (COE) 
is the space of the A-by-A symmetric unitary matrices S with the probability measure that 
satisfies the following invariance condition: 

(d H S) = {d H V T SV) 

for arbitrary unitary V. 

The circular symplectic ensemble (CSE) is the set of all A-by-A self-dual unitary matrices 
S with real quaternionic entries and with the probability measure that satisfies the following 
invariance condition: 

(d H s) = (d H J^V T J 2N Sv) , 

where S = <p (S) is the 2A-by-2A complex matrix representing S, V is an arbitrary unitary 
2A-by-2A matrix, and J2N is as defined above (a block-diagonal matrix with 2-by-2 blocks 
if (j) on the main diagonal). For more details see Dyson's papers or Chapter 2 in Forrester's 

book ma. 
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The eigenvalues of the matrices from these ensembles are located on the unit circle and 
can be identified with angles 8k on the unit circle. The density of the eigenvalue distribution 
for these ensembles is given by 



, - TT < Ol < TT, 



en, n p-^ 

l<j<k<N 

with p = 1, 2, and 4 for COE, CUE, and CSE, respectively. (See Proposition 2.5 in iflOln . 

The eigenvalues form random point fields on the unit circle, which we denote ', 
X£ UE , and X^ SE . For (3 = 2, there is a well known determinantal formula for the correlation 
functions of the random point field, and therefore X^ UE is determinantal. For (3 = 1 and 4, 
Dyson discovered formulas for the correlation functions in terms of quaternion determinants. 
Therefore, fields X^ OE and X^ SE are quaternion determinantal, and we can seek to apply 
the results which we found above. 

In this paper we will be concerned only with the symplectic ensemble. 

Let iV be any positive integer. Define 

1 a\ 1 sin(iV0/2) 



2vr ^ it ^ v ' 2ir sin (0/2) 

p p>0 v ' ' 

(Here the first summation is over (-2N + 1) /2, (-2N + 3)/2, (2N - 1) /2, and the 
second summation is over 1/2, ... , (2N — 1) /2.) Note that S2N (0) is even in 8. 
We write 

Ds 2N (8) := ^s 2N (8) = ± J>e^ e = ~\ 5>sin(p0) , 

p p>0 

and 

Is 2N (0) := / s 2N {8') d8', 
Jo 

so that 

IS2N (0) = ^-Y^P' 1 ^ 6 = 1 E - Sin ^ 
ZTTl TT p 

p p>0 

The functions Ds2N, and Is2N are odd in 8. 

Define the quaternion function a jv4 (0) by its matrix representation: 

(8) = -( S2N ^ Ds ^(0)\j_^( cos (p8) -p sin (p8)\ 
aN4[ ' 2\Is2n(0) s 2N (e) J 2tt ^[p- 1 sin (p9) cos ( P 8) )' 

In terms of quaternions and slightly abusing notation, the kernel can be written as follows: 



a N4 (8) = - ^s 2 N ~ - (Is2N + Ds 2 n) (Is2N - Ds 2 n) j 

N-l/2 

— \ (cos p8 + a p sin p8), 
2tt — ' 



p=l/2 
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where 

1 



ap= 2^ p2 - 1 )* i+ ^ 2 + 1 ) j ]- 
It is easy to check that cfy = — 1. 

Dyson proved the following result (Compare Theorem 3 in (H): The random field X^ SE 
is quaternion determinantal with the kernel <t/v4 (0 — 0')- 

The kernel a^i is a projection on a subspace of a finite-dimensional linear space L, where 

\N + 3 2N — 1 



L = span j 




-2N + 1 
2 








= span j 


— ^ cos p 


0, — = sin 









2 ' ' 2 

_ 1 3 2N - I 
P ~ 2' 2'"'' 2 
Note that dim Q L = 2N. 

If we calculate the action of the kernel a^A (0) in the basis 7r -1 / 2 cospO, 7r -1 / 2 sinpO, 
then we find the following matrix representation of the operator with kernel <7jv4 : 

..ef 1 a »-w)\. 

\ -OjV-l/2 1 / J 

(22) 

That is, if is a 2iV-by-2iV block-diagonal matrix that has blocks 





on its main diagonal. 

It follows that it can be written as 

v 

where 



v p = (°> • • • '°' i' ?' ' • • • >°) 



l 

V2 

and the entries 1 and a p are at places 2p and 2p + 1, respectively. 

Hence the kernel a^i (0) has a quasi-real diagonal form with eigenvalues A p = 1. It is 
easy to check that all iV-by-iV minors of the matrix K in ((22b have the determinant equal to 
either or 1, which implies that the kernel is completely positive. 

It follows from Theorem 1 1.41 that kernel a^b (#) defines a valid determinantal field. Of 
course, we already know this from Dyson's result. However, it is noteworthy that we proved 
our result by a different method, without relating the kernel to the eigenvalue distribution of a 
random matrix ensemble. 

Clearly the restriction of this field to an interval I = (a, b) is also quaternion determinantal 
with the kernel &iNb = 1/ (0) 0"JV4 (0 — 0') 1/ (#')• Numerical evaluations suggest that this 
restricted kernel is also quasi-real with positive eigenvalues for arbitrary N and I. However, 
the proof of this claim is elusive. The author proved it only for N = 2. 
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Appendix A. Proof of Theorem I3.1I 

Since the identity is algebraic, it is enough to show that it holds for matrices with real 
quaternion entries. We will prove this by showing that the corresponding result holds if we 
write the quaternion determinants in terms of Pfaffians. Namely, let C = X\ + X%i + X3J + 
X^k and define complex matrices A = X\ + X%i and B = X% + X&i. Then, by using ([8]) we 
obtain that 



-B*A + {B*A) t (A*A) l + B*B 



Det M (C*C) = -Pf 



-A*A-{B*B) t A*B-(A*B st 



The blocks -B*A+(B*Af and A*B-(A*B) t are antisymmetric, and ((A*A) t + B*B) = 
—A* A — (B*B) , so the block matrix is antisymmetric as well. 
What we need to prove is that 

/ -B*A + {B*Af (A*Af + B*B \ _ / -B I *A I + (B^A 1 )* (A 1 * A 1 )* + B I *B I 

\ -A* A - (B*B) t A*B - (A*B) t / j \ -A 1 * A 1 - {B u B ] ) t A u B l - (A^B 1 )* 

where the summation is over all ordered m-tuples I = (ii < . . . < i m ) , with € {1, . . . , n}, 
and A 1 , B 1 are the matrices that are obtained from matrices A and B, respectively, by taking 
the rows with indices in I. 

In order to prove this, we recall that if R is a 2m-by-2m antisymmetric matrix, then the 
pfaffian of R is defined as follows: 



2 m m! 



Pf ( i? )=« E ^n(a)n^(2,-iM 2i ). (23) 
Next, we note that if 

R = 



-B*A+{B*A) 1 (A*A) l + B*B 
-A*A-{B*B) t A*B-(A*bY 



then there is a formula for Rij in terms of elements of A and B. This formula depends on 
whether i and j are greater or less than m. For example if i and j are both < m, then 

n n 
a=l a=l 

If we substitute this in formula d23l . and expand, then we get 

Pi(R) = ^- i E *<7n(*) J] *a(^(l)^(2))...^(a(2m-l),a(2m)), 

o"GS2m (a,fe,...,z) 

(24) 

where the summation is over all m-tuples (a, 6, ... , z) with each letter taking a value in 
{l,...,n}. 

A similar formula holds for the pfaffian of R 1 , where 

R i _( —B 1 * A 1 + (B^A 1 ) 1 (A 1 * A 1 ) 1 + B^B 1 \ 



-A 1 * A 1 - (B^B 1 ) A T *B J - (A^B 1 ) 
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Namely, 

Pf ^ = vkri S E *« (°- C 1 ) ' CT ( 2 )) • • • ( ff ( 2m 1) ) ^ (2m)) . 

' cre5 2m (a,b,...,z)G/ m 

(25) 

The difference with the previous formula is that the elements of m-tuples (a, 6, ... , z) are 
now restricted to take values among indices in m-tuple I. 

Let us for shortness write *ff a> ... )Z (o") for the product \P a (<r (1) , a (2)) . . . ^ z (a (2m — 1) , a (2m)) 
If all elements of (a, 6, ... , z) are different, then the term ^ a ,...,z (c) occurs once in ex- 
pansion (l24l and once in the sum of expansions (|25| ). 

E pf M- 

(In this sum, it occurs in the expansion of that Pf (R 1 ) , for which / is the ordered version of 
the m-tuple (a, b, ... , z) .) 

If some of the elements of (a, b, ... , z) coincide, then the situation is different. The term 
^a,...,z (c) occurs once in expansion (1241 but it can occur more than once in the sum 

E pf M- 

For example, if all elements of the m-tuples are the same, a = b = . . . = z, then this term 
will appear in the expansion of each Pf (R 1 ) , whose index I contain a. 

Clearly, in order to prove that Pf (R) = Y^i Pf > ft * s enough to prove that the sum 
of all these terms is zero. That is, it is enough to show that for a fixed m-tuple (a, b, ... , z) 
with at least two elements that are equal, the sum 

sgn{a)^ a ^ z (a) 

is zero. 

Without loss of generality we can assume that a = b. We have to consider several cases of 
a, which are summarized in the following table 



a (I) 


a (2) 


a (3) 


a (4) 


sWi [i] 


SWrt [2] 




* 


* 


SWn [3] 


* 




* 


S 2m [4] 


* 


* 




52m [5] 






* 


S2m [6] 


* 






S 2m [7] 




* 




SWn [8] 


* 


* 


* 



The star means that the corresponding a (i) is greater than m. For example, S^m [1] denote 
the set of all permutations from S2 m th at satisfy the condition that all of a (1) , a (2) , a (3) , a (4) 
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are smaller than or equal to m. For permutations in this set, 
tf a (a (1) , a (2)) * Q [a (3) , a (4)) = {-B a4 A aJ + B ad A a4 ) (-B a , k A a j + B a jA a)k ) . 

S2m [2] denote the set of all permutations from S^m that sastify the condition that a (1) , a (2) 
are smaller than or equal to m and <r (3) , a (4) are greater than m, and so on. 

Let us define t\ [a] , as a permutation that coincides with a on all indices except 2 and 
4, for which it is defined by equalities t\ [a] (2) = a (4) , and t\ [a] (4) = <r (2) . Similarly, 
T2 [cr] is defined as a permutation which acts on everything except 2 and 3 as a, and on these 
indices it is defined by T2 [a] (2) = cr (3) , T2 [cr] (3) = a (2) . Finally we define T3 [<r] as a 
permuation that coincides with cr on all indices except 2, 3, and 4, where it is defined by the 
rules: r 3 [a] (2) = a (4) , r 3 [a] (3) = a (2) , r 3 [a] (4) = a (3) . Note that sgn (n [a]) = 
sgn (Y2 [cr]) = — sgn (a) , and sgn (r 3 [a]) = sgn (cr) . Observe that for an arbitrary function 
/, 

E (/M + /(nM) + /faM)) = 3 £ /(a). 

It is easy to check that the identity 

= * a (a (1) , a (2)) (a (3) , a (4)) - * a ( Tl [a] (1) , n [<r] (2)) * a (n [a] (3) , n [a] (4)) 
-* a (r 2 [a] (1) , r 2 [cr] (2)) * a (r 2 [a] (3) , r 2 [cr] (4)) 

holds for permutations in S2 m [1]» an d this implies that 

E ssn (a) (cr) = 0. (26) 

Next, 

E (/W + /(nH) + /(7s[a]))= e /( ct ) + 2 E /(*)■ 

and it is easy to check the identity 

= * a (a (1) , a (2)) \|> a (a (3) , a (4)) - * a ( n [a] (1) , n [a] (2)) * a (n [a] (3) , n [a] (4)) 
+* a (r 3 [a] (1) , r 3 [cr] (2)) * a (r 3 [a] (3) , r 3 [cr] (4)) . (27) 

for cr 6 S^m [2] . Hence, identity (1271 implies that 

E +2 E I <28) 

v o-eS 2m [2] cres 2m [3] 

The other cases are similar. We use 



E (f(a) + f(T 1 [a]) + f(r 3 [a])) = 2 E /(*)+ E / 

o-eS 2m [4] o-G5 2m [4] rtS 2m [2] 
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in order to conclude that 

2 + E J s 5 re (°") ^a,..., 2 (o") = (29) 



By adding ([28j) and ([291), we obtain 

E + E + E )*^W*wW = 0- (30) 

v °-SS 2m [2] <^6S'2 m [3] aG5 2m [4]/ 

Next, the identity 

£ (/(a)+/(r 1 H) + /(r 3 H)) = 2 £ /(a) + £ /(a) 

0-£S2m[5] Cre52m[5] 0"G<S2m [6] 

implies that 

2 ^ + E ) ^(ff) («r) = 0, (31) 

0"GS2m[5] creS 2m [6]y 

and the identity 

£ (/(«T) + /(7iM) + /(7s[a]))= £ /(a) + 2 £ /(a) 

ceS2m[6] crGS2m[6] 0"6<S2m [7] 

implies that 

^ +2 J] Un(<r)* w (a) = 0. (32) 

\ k cr65 2m [6] a-GS 2m [7]/ 

By adding (f3TT > and d32l >, we obtain: 

E + E + E ) sgn(a)^ a ,..., z (a) = 0. (33) 

^f£S2m[5] crGS2m[6] cr&S2m[t\/ 

Finally, 

E (/W + /(nM) + /(^H))=3 £ f(*). 

ves 2 m[8] <reS2m[8] 
Identity (|27T ) still holds in this case, and therefore, 

^2 sgn(a)* a _ z (a) = 0. (34) 

o-GS 2m [8] 

Clearly, the sets S2 m [k] , are disjoint and their union over k = 1, . . . , 8 consists of all 
permutations from S^m for which a (1) < m. Hence, identities (l26l ). (|30T >. (l33l . and (l34l 
imply that 

^ s#n (a) ^ a> ... >z (a) = 0. 

CGSWl 

er(l)<ra 

The case a (1) > m can be handled similarly, and we obtain 

^ sgn {a) * 0) ..., z (<r) = 0. 

0"GS2m 
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This holds provided the first two indices in (a, ... , z) are equal. It is clear that this identity 
also holds if any two indices in (a, . . . , z) are equal. 

As observed earlier, this implies that we can remove all terms *& a ,...,z (c) j f° r which two 
indices in (a, ... , z) coincide, from the expansions of both Pf (R) and the sum J2i Pf • 
(See expansions (l24l) and (|25T>.) Then it is clear that the remaining terms in the expansions are 
the same. Hence 

Pf(i?) = ^Pf (R 1 ), 
I 

and this implies the statement of the theorem. □ 
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